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ABSTRACT 

Most of the poorly known parameters of the Standard Model cannot be 
determined without reliable calculations in nonperturbative QCD. Lattice 
gauge theory provides a first-principles definition of the required functional 
integrals, and hence offers ways of performing these calculations. This paper 
reviews the progress in computing hadron spectra and electroweak matrix 
elements needed to determine as, the quark masses, and the Cabibbo- 
Kobayashi-Maskawa matrix. 
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1 Introduction 



Many contemporary reviews of elementary particle physics start by celebrating 
(or lamenting!) the success of the Standard Model. Indeed, with some nineteen^ 
parameters the SU(3)xSU(l)xU(l) gauge theory explains an enormous array of 
experiments. Even a terse compendium^ of the experiments is more than big enough 
to fill a phone book. A glance at Table || shows, however, that roughly half of the 
parameters are not so well determined. To test the Standard Model stringently, 
and thus to gain an inkling of what lies beyond, we must learn the values of these 
parameters more precisely. 

Except for the mass of the Higgs boson (or any other undiscovered remnant of 
electroweak symmetry breaking), the poorly known parameters all involve quarks. 
Other than topji which decays too quickly for confinement to play a role, the masses 
of the quarks are a bit better than wild guesses. The information on the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix is spotty, especially when one 
relaxes the assumption of three-generation unitarity, as shown in Table |2|. They are 
poorly determined simply because experiments measure properties not of quarks, 
but of the hadrons inside which they are confined. Of course, everyone knows what 
to do: calculate with QCD, the part of the Standard Model that describes the strong 
interactions. But then, the strong coupling is known only at the 5% level; not bad, 
but nothing like the fine structure or Fermi constants. Moreover, the binding of 
quarks into hadrons is nonperturbative — the calculations cannot be done on the 
back of an envelope. 

The most systematic technique for understanding nonperturbative QCD is lat- 
tice gauge theory. The lattice provides quantum field theory with a consistent and 
mathematically well-defined ultraviolet regulator. At fixed lattice spacing, the quan- 
tities of interest are straightforward (combinations of) functional integrals. These 
integrals can be approximated by a variety of techniques borrowed from statistical 
mechanics. Especially promising is a numerical technique, the Monte Carlo method 
with importance sampling, which has become so pre-eminent that the young and 
uninitiated probably haven't heard of any other. 

Results from lattice-QCD Monte Carlo calculations have begun to influence Ta- 
ble ||. The world average for the SU(3) gauge coupling as includes results from 
lattice calculations of the quarkonium spectrum,B i- H and at the time of this con- 
ference an even more precise result had appeared^ The same calculations are also 
providing some of the best information on the charmi and bottomi masses. This is 
an auspicious beginning. Over the next several years the lattice QCD calculations 
will mature. They will help to determine the other unknowns — light quark masses 
and the CKM matrix. The third column of Tables ^ and ^ lists relevant quantities 
or processes, and the rest of this talk explains how the program fits together. 

*I assume the neutrinos are massless, but count the "vacuum angle" of QCD. 
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Table 1: Parameters of the standard model and lattice calculations that will help de- 
termine them. Numerical values taken from the 1994 Review of Particle Properties,!! 
except rrit (Ref. §), sin 5, and 6'qcd- The strong coupling as refers to the MS scheme 
at Mz- Adapted from Ref. ^. 



parameter value or range related lattice calculations 
gauge couplings 

Qem 1/137.036 

W^Gf 1.166 GeV-2 

as 0.116 ±0.005 Amip_is; scaling 

electroweak masses 

mz 91.19 GeV 

mn > 58 GeV 
lepton masses 

me 0.51100 MeV 

ruf, 105.66 MeV 

m,r 1777 MeV 
quark masses 

vfiu 2-8 MeV m^, m\- 

md 5-15 MeV m^, rnj^ 

rus 100-300 MeV mj^ 

rric 1.0-1.6 GeV rnj/^ 

nib 4.1-4.5 GeV my 

mt 174 ± Wt\l GeV 
CKM matrix 

si2 0.218-0.224 K irev 

S23 0.032-0.048 B D*lu 

si3 0.002-0.005 B ttIu 

sin 5 /O Bk,Bb,Bb, 
QCD vacuum angle 

Oqcd < 10'^ dn 



Sect. ^ gives the non-expert some perspective on the conceptual and numerical 
strengths and weaknesses of lattice QCD. Sect. ^ reviews 1) the status of the light 
hadron spectrum and the propects for extracting niu, m^, and m^; and 2) results for 
the quarkonium spectrum, which yield as, rub, and rUc- Sect. |^ outlines lattice QCD 
calculations of electroweak, hadronic matrix elements that are needed to pin down 
the unitarity triangle of the CKM matrix. There are, of course, many other inter- 
esting applications to electroweak phenomenology; for more comprehensive reviews 
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Table 2: Ranges for CKM matrix elements \Vqr\ assuming unitarity but not three 
generations. Numerical values taken from the 1994 Review of Particle Properties. 
In three generations \Vud\ = si2, \Vcb\ = •S23, and \Vub\ = si3, to excellent approxi- 
mation. 



parameter value or range related lattice calculations 
0:974 

\Vus\ 0.218-0.224 

iKftI 0.002-0.005 

\Vcd\ 0.180-0.228 

\Vcs\ 0.800-0.975 

\Vcb\ 0.032-0.048 

\Vtd\ 0.0-0.13 

\Vts\ 0.0-0.56 

\Vtb\ 0.0-0.9995 



K — > TTeu 
B TTliy 
D ttIu 
D Klv 
B D*lu 



the reader can consult some of the papers listed in the bibliog raphyfl M 

2 Rudiments of Lattice Gauge Theory 

In quantum field theory physical measurements are related to vacuum expecta- 
tion values {O). Feynman's functional integral representation is 

(O) = lim limZ-^ f H d^.(x) TT #a(2;)#a(x) O e-^^^-^'^'^). (1) 

x,fi x,a 

The formula is easier to understand if read from right to left. S is the action — in 
our case the action of QCD, so it depends on the gluons and the quarks ip and 
anti-quarks ijj. O depends on the physics under investigation; the most useful kinds 
of O's are given below. The integration over all components and positions of the 
basic fields, with weight e*'^ would reproduce the familiar Schrodinger or Heisenberg 
formulations of quantum mechanics. The more convergent weight provides some 
benefits and imposes some restrictions — see below. Z is the same integral without 
O in it, so that (1) = 1. 

The limits are there to satisfy the mathematicians; without them the integrals are 
not well defined. These "cutoffs" also have a physical significance: we do not claim 
to understand physics either at distances smaller than a (the ultraviolet cutoff), 
or at distances larger than L (the infrared cutoff). The limit a ^ requires the 
renormalization group; it must be carried out holding L and physical, infrared scales 
fixed. In particular, the integration variables A^(x), il>{x), and ^(x) really represent 
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all degrees of freedom in a block of size a^. The limit "a — > 0" can be obtained not 
only literally, but also by improving the action of the blocked fields. 

These observations apply to any cutoff scheme for quantum field theory. A nice 
introduction to the renormalization-group aspects is a summer-school lecture by 
Lepagejni In lattice gauge theory a is nothing but the spacing between lattice sites. 
If there are on a side, L = Na. For given one can compute the integrals 
numerically. With the 10^-10^''-dimensional integrals that arise, the only viable 
technique is a statistical one: Monte Carlo with importance sampling. 

To compute masses the observable O = <I>(t)<I>^(0), where <I>(t) is an operator 
at time t with the flavor and angular-momentum quantum numbers of the state 
of interest. One can construct such operators using symmetry alone. The radial 
quantum number would require a solution of the theory, but that's what we're 
after. The "two-point function" 

i'^mHo)) = Y.\m\n)\\-^"\ (2) 

n 

where the sum is over the radial quantum number. The exponentials are a happy 
consequence of the weight in eq. It is advantageous because at long times t 
only the lowest-lying state survives. In a numerical calculation masses are obtained 
by fitting two-point functions, once single-exponential behavior is verified. Since 
$ is largely arbitrary, some artistry enters: if single-exponential behavior sets in 
sooner, the statistical quality of the mass estimate is better. 

To compute a matrix element of part of the electroweak Hamiltonian, TC, the 
observable O = $7r(i7r)'^(ih)^s('-') transition from hadron "i?" to hadron 

"tt." At long times th and — th the "three-point function" 

plus excited-state contributions. If, as in decays of hadrons to leptons, the hadronic 
final-state is the vacuum, a two-point function will do: 

(nmUo)) = E(0|W|i?„)e-™"*(i?„|<I>l,|0). (4) 

n 

The desired matrix elements {Tr\7i\B) and {0\7i\B) can be obtained from eq. (|3|) 
and because the masses and <I>-matrix elements are obtained from eq. (^). 

To obtain good results from eqs. (|2|)-@, it is crucial to devise nearly optimal op- 
erators in the two-point analysis. Consumers of numerical results from lattice QCD 
should be wary of results, still too prevalent in the literature, that are contaminated 
by unwanted states. 

In the numerical work that mostly concerns us here, the integrals are computed 
at a sequence of fixed a's and L's. One adopts a standard mass, say rup, and defines 



expt 

rup ^ 



(5) 
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to obtain the latice spacing in physical units, and other quantities are predicted via 

mB = . (6) 

a 

For continuum-hmit, infinite- volume results this is the same as extrapolating dimen- 
sionless ratios, e.g. 

^ = lim lim <'^B{a,L) 
nip L^ooa^o airipla, L) 

There is theoretical guidance for both limits. According to general properties 
of massive quantum field theories in finite boxes ,0 the infinite- volume limit is 
rapid for niT^L ^ 1, exponential for masses. In non-Abelian gauge theories the 
renormalization-group a ^ limit is controlled by asymptotic freedom. 

The main strength of lattice QCD is that it is QCD. It has 1 + Uf adjustable 
parameters, corresponding to the gauge coupling and the quark masses. From the 
renormalization group, the adjustment of the gauge coupling is equivalent to setting 
the lattice spacing in physical units, eq. (^). Once the parameters are determined 
by 1 + nf experimental inputs, QCD should predict all other strong- interaction 
phenomena. There is no need to introduce condensates, as in ITEP sum rules, 
or non-renormalizable couplings, as in chiral perturbation theory or heavy-quark 
effective theory. If theory and experiment disagree, it is a signal of new physics. 

There are some disadvantages. A practical, though not conceptual, problem is 
that large-scale computational work is more labor-intensive than traditional theo- 
retical physics. Careful work is needed to estimate the uncertainties reliably. The 
improvements in computer power and algorithms of recent years have helped practi- 
tioners understand their uncertainties better and better. As the consumers of their 
results become commensurately sophisticated, this trend will continue. After all, in 
the context of Table |l|, meaningful error bars are just as important as the central 
value. 

2.1 The quenched approximation 

The biggest disadvantage of most of the numerical results mentioned in this 
talk is something called the "quenched" approximation. A meson consists of a 
valence quark and anti-quark exchanging any number of gluons. The gluons can 
turn into virtual quark loops and back again. The latter process costs a factor of 
100-1000 in computer time, so many Monte Carlo programs just omit the virtual 
quark loops. To accentuate the positive — gluons and valence quarks are treated 
better than in non-QCD models of hadrons — the omission is sometimes called the 
valence approximation. To admit the negative, it is less often called the loopless 
approximation. But most often lattice mavens borrow a jargon from condensed- 
matter physics and call it the quenched approximation. 

If quenched QCD makes any sense, it is as a kind of model or effective theory. 
The parameters of quenched QCD can be tuned to reproduce physics at one scale. 
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But the /? function of quenched and genuine QCD differ, as one sees in perturbation 
theory, so one cannot expect agreement at ah scales. As with any model, only in 
special cases can one argue that these effects are correctable or negligible; these 
cases will be highlighted in the rest of the talk. 



3 From Hadron Spectra to the QCD Parameters 

3.1 Light hadrons and light-quark masses 

Over the past few years a group at IBM has carried out a systematic calculation 
of the light-hadron spectrum using the dedicated supercomputer GFll.lli They have 
numerical data for 5 different combinations of (a,L). At L ~ 2.3 fm there are 
three lattice spacings varying by a factor of ~ 2.5. At the coarsest lattice spacing 
(a~^ ~ 1.4 GeV) there are three volumes, up to almost 2.5 fm. A variety of quark 
masses are used, and the physical strange quark is reached by interpolation, whereas 
the light (up and down) through extrapolation. The mass dependence is assumed 
linear, as expected from weakly broken chiral symmetry, and the data substantiate 
the assumption. 

The units (i.e. lattice spacing) has been fixed with nip and the quark masses 
with and rrix- The final results, after extrapolation to the continuum limit and 
infinite volume, are shown in Fig. || for two vector mesons and six baryons. (The 
quark- mass interpolation could reach only the combination mE+m-s — m^.) Despite 
the quenched approximation the agreement with experiment is spectacular. 

Fig. H also includes results from the same investigation for decay constants. 
The agreement of f-n/rup and fK/inrip is not as good as for the masses. Because of 
the quenched approximation, this is not entirely unexpected. Recall the argument 
concerning distance scales and effective theories in sect. |2.1| . The binding mechanism 
responsible for the masses encompasses distances out to the typical hadronic radius. 
The decay constant, on the other hand, is proportional to the wavefunction at 
the origin and thus is more sensitive to shorter distance scales. One sees better 
agreement when forming the ratio fK/f-m which — recall eq. (^) and subsequent 
discussion — is like retuning to the shorter distance. 

One would like to use the hadron masses to extract the quark masses. Because 
of confinement, the quark mass more like a renormalized coupling than the classical 
concept of mass. Calculations like the one described above yield immediately the 
bare mass of the lattice theory. More useful to others would be the MS scheme 
of dimensional regularization. A one-loop perturbative calculation can be used to 
convert from one scheme to another Jll' O i 

For the light quarks it is convenient to discuss the combinations rh = ^(mrf+mu), 
^"^dii ~ "^d ~ "^u' nis. Ratios of the light-quark masses are currently best 
estimated using chiral perturbation theory.0 To set the overall scale requires a 
dynamical calculation in QCD. In lattice QCD, rh and can be extracted from 
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Figure 1: The spectrum and decay constants of the hght hadrons. Error bars 
are from lattice calculations in the quenched approximation, and • denotes 
experiment. 



the variation in the square of the pseudoscalar mass between and m]^. The 
most difficult quark-mass combination is Am^^, which causes the isospin-violating 
end of the splittings in hadron multiplets. Since chiral perturbation theory provides 
a formula for Am^^/m^ with only second-order corrections, it is likely that the best 
determination of Am^^ will come from combining the formula with a lattice QCD 
result for m<j. 

Using the compilation of quenched and unquenched results of Ukawajll Mack- 
enzieii has estimated mj^(l GeV) ~ 2.3 MeV and rn^y^{l GeV) ~ 65 MeV. The 
symbol ~ stresses the lack of error bar. This is outside the ranges of 3.5-11.5 MeV 
and 100-300 MeV indicated in Table 111. A more recent analysis of the strange quark 
finds my^^{2 GeV) = 127 ± 18 Mev|3 in the lower part of the range in Table [l|. 
None of these results should be taken seriously until a more complete error analysis 
exists, but it is intriguing that the conventional estimates might be too high. 
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3.2 Quarkonia, as, and heavy-quark masses 

Quarkonia are bound states of a heavy quark and heavy anti-quark. Three 
famihes of states exist, charmonium (r/c, J/ip, etc), bottomonium (r^^, 'Y^, etc), and 
the as yet unobserved Be {be and be bound states). Compared to hght hadrons, 
these systems are simple. The quarks are nonrelativistic, and potential models give 
an excellent empirical description. But a fundamental treatment of these systems 
requires nonperturbative QCD, i.e. lattice QCDJ^ Potential models can be exploited, 
however, to estimate lattice artifacts, and in the quenched approximation they can 
be used to make corrections. Many states have been observed in the lab, providing 
cross-checks of the methodology of uncertainty estimation. 

Once the checks are satisfactory, one can use the spectra to determine as, Trie, 
and mf). One can also have some confidence in further applications, such as the 
phenomenology of D and B mesons discussed in sect. 

For charm, and especially for bottom, the quark mass is close to the ultraviolet 
cutoff, 1/a or vr/a, of present-day numerical calculations. Originally lattice gauge 
theory was formulated assuming m^a <C 1, so quarks m^a ~ 1 require some reassess- 
ment. There are four ways to react. The patient, stolid way is to wait ten years, 
until computers are powerful enough to reach a cutoff of 20 GeV — not very inspir- 
ing. The naive way is to extrapolate from smaller masses, assuming the m^a <^ 1 
interpretation of the lattice theory is adequate; history shows that naive extrapola- 
tions can lead to naive and, thus, unacceptable error estimates. The insightful way 
is to formulate an effective theory for heavy quarks with a lattice cutoff;^' H this 
is the computationally most efficient approach, and when the effectiveness of the 
heavy-quark expansion is a priori clear, it is the method of choice. The compulsive 
way to examine a wide class of lattice theories without assuming either m^a <C 1 or 
rUg ^> (AqcD)0~^); by imposing physical normalization conditions on masses and 
matrix elements, it is possible to interpret the correlation functions at any value 
of mga.il The underlying reason is that the lattice theory is completely compatible 
with the heavy-quark limit, so the mass-dependent interpretation connects smoothly 
onto both the insightful method for m^a ^ 1 and the standard method for niga <C 1. 

Fig. ^ shows the charmonium spectrum, on a scale appropriate to the spin- 
averaged spectrum. Light quark loops are quenched in these calculations.il' il The 
agreement with experimental measurements is impressive, but Fig. |2| barely displays 
the attainable precision. Fig. ^ shows the fine and hyperfine structure of the P states, 
now for bottomonium F^l (The ^Pi state hi, has not been observed in the lab; the 
he has been seen.) The authors of Ref. also have results with the virtual quark 
loops from two light quarks, i.e. up and down are no longer quenched, but strange 
still is. The agreement is comparable.l^ 

To obtain these results only two parameters have been adjusted. The standard 

^The utility of quarkonia for testing the methodology of lattice gauge theory and the theory of 
QCD has been stressed over and over by Peter Lepage. 
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Figure 2: A comparison of the charmonium spectrum as calculated in lattice QCD, 
using two different methods. Ref. |2^: o, Ref. p^: □. From Ref. & 



mass in eq. (^) is Amip_is, the spin-averaged splitting of the IP and IS states, which 
is insensitive to the quark mass. By the renormalization group, this is equivalent to 
eliminating the bare gauge coupling, or to determining Aqcd- The bare quark mass 
is adjusted to obtain the spin average of the IS states that is measured in the lab. 
Otherwise figs. § and |3| represent predictions of quenched QCD. 

The success of these calculations permits one to extract the basic parameters, 
as and m^. There are four steps: 

1. Compute the charm- and bottomonium spectra with ?t-/,mc = 0, 2 or 3 flavors 
of virtual quark loops. (nj^MC = corresponds to the quenched approxima- 
tion; nj MC = 2 quenches just the strange quark; nj^MC = 3 would be the real 
world.) 

2. With perturbation theory, convert the bare lattice coupling Qq"'^'^'~^^ to the 
quark-potential (V) or MS scheme; convert the bare lattice mass 

to the pole or MS scheme. The natural scale for this conversion is near (but 
not quite§i) vr/a. 

3. Unless ?i-/,mc = 3, correct for the quenched approximation. 
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Figure 3: Lattice QCD results for the spin-dependent splittings of the lowest-lying 
P states in bottomonium. The dashed lines are the experimental values, where 
available. Energies are measured relative to the spin average of the x states. From 
Ref. E5|. 



4. Eliminate a from aj^(7r/a) and amj^(7r/a) using 

_ aAmip-is 
" ~ 460 MeV ' ^ ^ 

where the numerator is the IP-IS splitting in lattice units. 

Steps H and ^ are explained above. Step |2| requires one-loop perturbation theory, 
suitably optimized.!! Step ^ is crucial, because without it the results have no busi- 
ness in Table |^. 

Consider first as, and recall the idea of treating the quenched approximation 
as an effective theory. One sees that the couplings are implicitly matched at some 
scale qg characteristic of quarkonia. So the matching hypothesis, supported by 
figs, m and ^, asserts 

«^"'^('?q) = 4'^(^q)- (9) 
Potential models teh us that 200 < Qc < 800 MeV and 200 < % < 1400 MeV. Step 
3 yields a^^'^'~^^(7r/a), so one can use the two- loop perturbative renormalization 
group to run from vr/a to qq. The perturbative running is an overestimate if qq is 
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taken at the lower end of these ranges .[] This argument was used for the original 
lattice determinations of the strong coupling, and its reliability was confirmed in 
^/,MC = 2 calculations .§1 

Currently the most accurate result is from Ref. 0, 

aPi8.2 GeV) = 0.196 ± 0.003, (10) 

based on re/,MC = and n/,MC = 2 results, with an extrapolation in nj. The V 
scheme is preferred for the matching argumerit, not only for physical reasons, but 
also because of its empirical scaling behavior.c^ The scaling behavior implies that 
one can run with the two-loop renormalization group to high scales and convert to 
other schemes. For comparison to other determinations, eq. (|l^) corresponds to 

ay^{Mz) = 0.115 ± 0.002. (11) 

The quoted uncertainty is smaller than that reported from any other method. The 
largest contributor is the quenched correction; the second largest is the perturbative 
conversion — > V ^ MS. 

To determine the quark mass the one applies the same renormalization-group 
argument. But, quark masses don't run below threshold!! Hence, for heavy quarks^ 
"^Q^'^^Vg) = "iQ''"^^('?Q) = m§\QQ)^^'i\rnQ). The only corrections are 
perturbative, from lattice conventions to MS or physical conventions. Using the 
convention of the perturbative "pole mass" 

nic = 1.5 lb 0.2 MeV (Ref. 8, preliminary), 

mb = 5.0±0.2 MeV (Ref. |). 

At the nonperturbative level confinement wipes out the pole, so the name "pole 
mass" should not be taken too literally. The perturbative pole mass is like a running 
mass, except that it has a fixed scale built into the definition. It is useful, because 
it is thoug htS to correspond to the mass of phenomenological models that do not 
probe energies less than Aqcd; such as potential models. In other contexts — such 
as the study of Yukawa couplings in unification scenarios — the MS convention may 
be more appropriate. Eq. (p^) corresponds to Ms("^^) ~ ±0.1 GeV. 

4 From Matrix Elements to the CKM Matrix 

Electroweak decays of fiavored hadrons follow the schematic formula 



experimental 
measurement 



known 
factors 



QCD \ / CKM 
factor / I factor 



(13) 



*For light hadrons, giight ~ Aqcd, so there would be no perturbative control whatsoever. 
^For light quarks the threshold is deep in brown muck, and all bets are off. 
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North American, Japanese, and European taxpayers provide us with lots of money 
for the relevant experiments, because they want to know the CKM factors. But 
unless we calculate the inherently nonperturbative QCD factor, they will be sorely 
disappointed. 

It is convenient to start with the assumption of three-generation unitarity. Then 

VudV:^ + + VtdVtl = 0, (14) 

an equation that prescribes a triangle in the complex plane. Dividing by VcdVcb 
writing VudVub/^cdVd) = P + one sees that unitarity predicts 



VcdV 



ir]. (15) 



cb 



The notationpO| (/j, fj) is to distinguish these parameters from the standard Wolfen- 
stein parameters (/O, r/) = {p,fj)/\Vud\- The standard CKM phase in Table || is 
6 = tan~^ {rj / p) = tan~^(ry//)). 

One would like to determine (p, ff) through as many physical processes as pos- 
sible. For example, the strength of CP violation in 6-flavored hadrons is related to 
the angles of the triangle, hence the high interest in B factories. If all experiments 
agree, the test verifies the CKM explanation of CP violation; if not, the discrepency 
would have to be explained by physics beyond the Standard Model. Meanwhile, lat- 
tice QCD is useful for measuring the sides of the unitarity triangle. Depending on 
the shape of the triangle, the precision may be good enough to predict the angles 
before the B factories have been built. 

Let us consider the CKM matrix elements in eq. (^) . Assuming three-generation 
unitarity, 1 — \Vtb\ is too small to worry about, and \Vcd\ — \Vus\ is also very small. 
In principle, however, \Vus\ and \Vcd\ can be determined from lattice QCD and 
measurements of the semi-leptonic decays K iiev and D ttIv, respectively. 
The technique is the same as for \Vub\ from B — > nlv^ discussed below. It is unlikely 
that lattice QCD can, or will need to, improve on \Vud\ = 0.9744 it 0.0010 during 
the period relevant to this discussion. 



The most poorly known elements of eq. (14) are \Vcb\^ \Vub\i and \ Vtd\- In prin- 
ciple, the first two can be determined from leptonic decays Bq — > tv, q = u, b, but 
the experimental prospects are bleak. The semi-leptonic decay is more promising. 
Near g^ax = i^B — rriD*)'^ the differential decay rate for B — > D*lv is 



647r^mB 



imax 



Ax{q')\'\Vcb\\ (16) 



where is the invariant mass-squared of the lepton system. One must carry out 
a nonperturbative QCD calculation to obtain the form factor A\((j^\ By heavy- 
quark symmetry, however, ^i((i'max) obeys a normalization condition, up to l/m|). 
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Figure 4: The Isgur-Wise function S,{uj) (essentially the form factor Ai of the text) 
from lattice QCD and CLEO. The kinematic variable uj = Va • Vb = 1 — ((?max ~ 
q'^)/2mB'm£)* . From Ref. 36 



corrections^ (estimated to be small) and known radiative corrections. Other form 
factors, which are phase-space suppressed near (7^ax' ^^^o related by heavy-quark 
symmetry to ^i(g^). Hence, eq. (16) provides an essentially model- independent!! 
way to determine \Vcb\- 

The difficulty with the model-independent analysis is that the decay rate vanishes 
at (7max- To aid experimentalists' extrapolation to that point, several group 

11,0 i 

have used quenched lattice QCD to compute the slope of ^i. A typical analysis is 
to fit the slope to lattice-QCD numerical data, and then fit the normalization to 
CLEO's experimental data, as shown in Fig. ^. For example, Simone of the UKQCD 
Collaboration findsil 

The first error is experimental; the second is from the lattice-QCD slope. Unfor- 
tunately, it is not clear how to correct for the quenched approximation, and the 
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associated uncertainty has not been estimated. Moreover, consistency checks of 
varying lattice spacing, volume, etc, are still in progress. Nevertheless, the overall 
consistency with experiment, shown in Fig. ^ is encouraging. 

\Vub\ can be obtained from the semi-leptonic decays B plv and B — > ttIv. 



Expanding in near q^^-^ = {mB 
reads 



the differential decay rate for B — > t^Iv 



247r3 KmB/ 



2 

max 



\UiQ 



^ ) 

max/ 



|K6|'(l + 0(gLx-9' 



(18) 



where f+{q^) is the form factor that must be calculated in lattice QCD. Now, how- 
ever, heavy-quark symmetry does not restrict /+(gmax)) so a calculation is needed to 
make any progress. These calculations are underway at Fermilab, and, presumably, 
many other places. 

The third row of the CKM matrix can be probed via the box diagrams responsible 
for neutral meson mixing. The CP = + admixture of the is parameterized by 
|eii-| = 2.26 X 10^3. The Standard Model predicts 



IGvr^mxAmx 



\m\f'^BK\Vudyu. 



f?(|VcfeP(i - p)ytmh{yt) + ydmhivt) - m) 



(19) 



where yq = m'^/m'^r. This formula assumes three-generation unitarity and neglects 
the deviation of \Vcs\ and \Vtb\ from unity. The r/j and /i_multiplying the CKM 
factors arise from box diagrams and their QCD corrections .eZI The nonperturbative 
QCD factor is |m 

operator. The best result for Bk is^ 



j^ff^Bxi which is the K-K transition matrix element of a A5 = 2 



S/^(NDR, 2 GeV) = 0.616 ± 0.020 ± 0.014 ± 0.009 ± 0.004 ± 0.002 

± (few %) ± 3%. 



(20) 



The many error bars are exihibited to show how mature the uncertainty analysis 
has become. The first is statistical and the others are systematic. The "few %" are 
for the quenched approximation. This estimate comes from repeating some of the 
numerical computations for full qcd|o| though not enough to obtain the other error 
bars, and from an analysis of chiral logarithms.0 The latter study is reassuring only 
for degenerate quarks, so the calculations are done with both quarks at im^. The 
3% uncertainty is an estimate of 0(m<j — rrid) contributions. Combining the errors 
and converting to the renormalization-group invariant that appears in eq. (|l9|), one 
findsH 

^K = (aj^(/i))-6/255^(NDR, /x) 



:0.825±0.027(stat.) ±0.023(syst.) ± (few 
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This result places a high standard on calculations of B^, whether by lattice QCD 
or any other method. Would-be competitors must not only reach 10% uncertainties, 
they must do so with an error analysis as thorough and forthright as Ref. 

Mixing in the B^-B^ system is also sensitive to Vtd- In the Standard Model the 
mass splitting is given by 



Xd 



'mlJkB^lVCM', (22) 



3 i>d •' i'd i>d 



The same formula holds for the Bg-Bg system, but with the d quark replaced by s 
(i.e. Bd ^ Bg, Vtd ^ yts-) The nonperturbative QCD factor is ^m\^f'^^BBq, which 
is the Bq-Bq transition matrix element of a AS = 2 operator. The calculation of 
the decay constant /b has received a great deal of attention over the last several 
years,§i but the matrix element needed here, ^m%f%BB, has been mostly neglected. 



(There are some older, exploratory papers r^l) 

It is interesting to see how the lattice results influence the unitarity triangle. 
Fig. |5| shows constraints from eqs. (|l9|) , iVufc/Vcbl, and Xd/xg, taking for the masses 

"^c,MS = l-3±0.2, 
"^t,MS = l'^5±15, 

for the hadronic matrix elements 

Bi^ = 0.825 ±0.050, 
|/Bd//BJ= 0.90 ±0.05, (24) 
\BbJBbJ = 1.0±0.2, 

for "experimental''^] CKM results 

|Kb|= 0.040 ±0.005, 
IWKbN 0.08 ±0.02, ^ ^ 

and for neutral B mixing measurements 

Xrf = 0.72 ±0.08, 
x. = 15±5. 

Other inputs are as in Ref. I've made two wild guesses: \Bb^/BbJ and Xg- But 
note that I take the uncertainty estimate in Bk seriously; doubling it would not 
make much difference, in view of the uncertainties in nit and \Vcb\- Alas, these and 
the other uncertainties are too large to make Fig. ^ interesting. 

^Nonperturbative QCD is needed to extract these results! 
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Figure 5: Constraints on (p, f/) from {ekI (solid hyperbolae), l^fe/Fctl (dashed circles 
with origin (0,0)), and Xd/xs (dash-dotted circles with origin (1,0)), and contempo- 
rary uncertainties. 



What if lattice QCD calculation and the experiments improve? Consider for the 
masses 

'^c,MS = l-3±0-l' .27) 
"^t,MS = 1 ''S ± 5, 

for the hadronic matrix elements 

= 0.825 ± 0.027, 
|/b,//bJ= 0.90 ±0.02, (28) 

in particular eliminating almost all the statistical error in 13 k', for "experimental" 
CKM results 

1^61 = 0.035 ±0.002, 
|Kfe/V'cf,| = 0.080 ±0.004 "low," (29) 
= 0.091 ±0.004 "high," 

and for neutral B mixing measurements 

Xrf = 0.72 ±0.04, 

X, = 18 ± 2. ^'^^^ 
Fig. ^ shows how this 5-10% level of precision improves the limits on {p,f]). The 
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Figure 6: Constraints on (p, 77) from {ekI (solid hyperbolae), low |14f,/V"cfe| (dashed 
circles with origin (0,0)) or high |Ktfe/K:f)| (dotted circles with origin (0,0)), and 
Xd/xs (dash-dotted circles with origin (1,0)), and improved (5-10%) uncertainties. 

wildest guess remains Xg, so ignore the dashed-dotted curves momentarily. The 
region allowed by the hyperbolic band from ek and the circular band from 
shrinks if \Vub\ is too small. The tension between these two constraints is partly a 
consequence of the low value of \Vcb\ suggested by eq. (|l7|). Increasing \Vcb\ brings 
the hyperbolic band down more rapidly than it shrinks the circular band. 

If the real- world values of \ Vcb\ and iKifc/Vcbl allow a sizable region, as for the dot- 
ted circles in Fig. neutral B mixing becomes crucial. The constraint becomes more 
restrictive as Xg increases. Unfortunately, the experimental measurement becomes 
more difficut as Xs increases. If it proves impossible to obtain useful information on 
Xg, one can return to eq (^), and focus on x^ alone. The lattice-QCD calculations 
of Im^/g^B will carry larger uncertainties, however, than the : Bg ratio. 

5 Conclusions 

This talk has examined several ways in which lattice QCD can aid the deter- 
mination of standard-model couplings. The quenched lattice calculations may be 
divided into several classes, according to the maturity of the error analysis and the 
presumed reliability of the quenched approximation. One class consists of the light- 
hadon and quarkonia spectra and the K-K mixing parameter Bk- For them the 
straightforward uncertainties (statistics, a, L, excited states, perturbation theory) 



17 



seem fairly estimated. The quenched approximation is another matter. In quarko- 
nia, one can correct for it with potential models, yielding determinations of as and 
the charm and bottom masses. The quenched error in Bk is also thought to be 
under control, and — taking the error bars at face value — Bk is no longer the limit- 
ing factor in the le^l constraint on the unitarity triangle. A second class consists 
of /b, the semi-leptonic form factors of K and D mesons (not discussed in this 
talk, but see Ref. 42), and the Isgur-Wise function. These quantities are essential 
for direct determinations of the first two rows of the CKM matrix. The quenched- 
approximation calculations are in good shape, but the the corrections to it cannot 
be simply estimated. A third class consists of the semi-leptonic decay B ttIu and 
neutral B mixing, for which only exploratory work has appeared. 

Nevertheless, all QCD quantities discussed here will follow a conceptually clear 
path to ever-more-precise results. The next ten years or so will almost certainly wit- 
ness computing and other technical improvements that will allow for wide-ranging 
calculations without the quenched approximation. By then the most efficient tech- 
niques for extracted the most relevant information will have been perfected. 
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